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Abstract 

We have previously shown that a very small amount of Lorentz invariance vi- 
olation (LIV), which suppresses photomeson interactions of ultrahigh energy 
cosmic rays (UHECRs) with cosmic background radiation (CBR) photons, 
can produce a spectrum of cosmic rays that is consistent with that currently 
observed by the Pierre Auger Observatory (PAO) and HiRes experiments. 
Here, we calculate the corresponding flux of high energy neutrinos generated 
by the propagation of UHECR protons through the CBR in the presence of 
LIV. We find that LIV produces a reduction in the flux of the highest energy 
neutrinos and a reduction in the energy of the peak of the neutrino energy 
flux spectrum, both depending on the strength of the LIV. Thus, observa- 
tions of the UHE neutrino spectrum provide a clear test for the existence 
and amount of LIV at the highest energies. We further discuss the ability of 
current and future proposed detectors make such observations. 
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1. Introduction 



Ultrahigh energy cosmic rays and neutrinos are of interest as possible 
probes of new physics [l|. In particular, some quantum gravity models pre- 
dict that Lorentz invariance may be weakly broken at the very high energies, 
leading to potentially observable consequences. The possibility of using ultra- 
high energy cosmic rays (UHECRs) to probe for a small violation of Lorentz 
invariance was suggested over a decade ago 0]. Indeed, a detailed analysis 
of the effects of LIV on the UHECR spectrum has yielded the tightest con- 
straint on LIV to date [jjj]. Shortly after the discovery of the CBR it was 
pointed out that photomeson interactions of UHECRs with photons of the 
cosmic background radiation (CBR) would result in a sharp steepening of 
their spectrum above E ~ 50 EeV now known as the "GZK effect" jj[5[. 
However, even a very small amount of LIV will kinematically inhibit some of 
these interactions. It has been previously shown that a possible signature of 
LIV in the UHECR spectrum would be a recovery of the cosmic ray spectrum 
at energies greater than ~ 200 EeV (0, @]). 

Given the current state of the UHECR observational data as reported by 
HiRes 0, 0] and Auger Q, it is possible to constrain LIV, but not to rule 
it out. However, unlike cosmic-ray baryons and photons, ultrahigh energy 
neutrinos do not suffer significant energy losses over cosmological distances. 
Studies of UHE neutrinos with both ground-based and space-based detectors 
could provide a new and less ambiguous test of LIV. This is because they 
are a guaranteed byproduct of the photomeson interactions of UHECRs with 



the CBR followed by subsequent pion decay [lOj, [ll|, |12J, |13 



In this paper, we further consider the observational implications of the 
effect of a very small amount of LIV, viz. the suppression of photomeson 
production on the resulting UHE neutrino spectrum. In our previous work 
pl, [sj we undertook a detailed calculation of the modification of the UHECR 
spectrum caused by LIV using the formalism of reference ^ and the kine- 
matical approach originally developed in reference 14j. We employ the same 
techniques used in references [i[ and to determine the resulting photome- 
son neutrino spectrum. We again consider here the case where the primary 
UHECRs are protons. 

ANITA II 15[ has placed an upper limit on the neutrino flux for energies 
greater than > 10 18 eV. IceCube is nearing completion and will be sensi- 
tive to neutrinos of energies up to ~ 10 17 eV. Proposed future ground-based 
and space-based neutrino detectors could be capable of detecting and study- 



2 



ing photomeson neutrinos. We will discuss the ability of such detectors to 
constrain LIV or observe its effect. 

2. LIV and the Spectrum of UHECRs 

We now extend the calculation of reference [sj to determine the photome- 
son neutrino fluxesQ A full description of the LIV formalism we use is given 
in references ji[ and [3[. We summarize the salient points here. The free 
particle Lagrangian is modified by the inclusion of a leading order perturba- 
tive, Lorentz violating term. This term leads to the modified free particle 
dispersion relations 

E 2 = p 2 + m 2 + 25p 2 . (1) 
These relations can be put in the standard form 

E 2 = p 2 c 2 MAV + m 2 c 4 MAV , (2) 

by shifting the renormalized mass by the small amount m — > m/(l + 25) and 
shifting the velocity by the amount 

c MAV = y/(l + 25) -1 + 5 (3) 

where cmav is identified as the maximum attainable velocity of the free 
particle in the reference frame of the CBR. Using this formalism, different 
particles can have different maximum attainable velocities (MAVs) that can 
all be different from 1 as well as different from one another. Hereafter, we 
denote the MAV of a particle of type % by q and the difference 

d-Cj EE Si, (4) 

These modified dispersion relations are then applied to the kinematical 
relations governing the dominant single meson photomeson interaction: 

p + 7-^V + vr. (5) 

From equations §B) and (jlj), a dispersion relation can be constructed for a 
particle a 

E 2 =p 2 + 25 aP 2 + m 2 (6) 



We use the usual convention c = 1 . 
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where 5 a is the difference between the MAV for the particle a and the speed 
of light in the low momentum limit, c = 1. 

In order to modify the effect of photomeson production on the UHECR 
spectrum above the GZK energy, 5 wp > as shown in reference (2)]. The 
condition for photomeson interactions to take place is 

5 np < 3.23 x l(T 24 (u;M)) 2 - (7) 

where 10 is the energy of the CBR photon and oo = kT CB R = 2.35 x 1(T 4 eV 
with T CB R = 2.725 ± 0.02 K @, Q. 

If LIV occurs and 5 np > 0, photomeson production can only take place 
for interactions of CBR photons with energies large enough to satisfy equa- 
tion (J7J). This condition implies that while photomeson interactions leading 
to GZK suppression can occur for "lower energy" UHE protons interacting 
with higher energy CBR photons on the Wien tail of the spectrum, other 
interactions involving higher energy protons and photons with smaller values 
of to will be forbidden j^] . Thus, the observed UHECR spectrum may exhibit 
the characteristics of GZK suppression near the normal GZK threshold, but 
the UHECR spectrum can "recover" at higher energies owing to the possibil- 
ity that photomeson interactions at higher proton energies may be forbidden 
[3.1, |6|, ll6|, ll7J. Even a small violation of Lorentz invariance changes the in- 



elasticity of the interaction, (i.e., the amount of energy transferred from the 
incident proton to the created pion). This is the key to understanding the 
effect of LIV on photomeson production. With an increase in proton en- 
ergy, the range of kinematically allowed angles of interaction between it and 
the photon becomes more restricted, thus reducing the phase space and, in 
turn, the total inelasticity. Figure [fl reproduced from reference jij , shows 
the calculated inelasticity modified by LIV for a value of 5 np = 3 x 10 -23 as 
a function of both CBR photon energy and incident proton energy. Other 
choices for 5 np yield similar plots but change the energy at which LIV effects 
become significant. The inelasticity precipitously drops above a certain en- 
ergy because the LIV term in the pion rest energy from equation ([6]) becomes 
comparable to m n . 

The proton energy loss rate by photomeson production is given by 



1 dE ujqc 



E dt 27T 2 7 2 n 3 c 3 



00 

J dt e a{e)K(e) ln[l - e - e/2 ^ } (8) 
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Figure 1: The calculated proton inelasticity modified by LIV for S^p = 3 x fO 23 as a 
function of CBR photon energy and proton energy (from reference |3J). 



where e is the photon energy in the center of mass system, K(e) is the 
modified inelasticity calculated from the kinematics, and cr(e) is the total 7- 
p cross section. The lower limit of the integration, 7], is the photon threshold 
energy for the interaction in the center of mass frame. 

As in reference ^| , we assume that the source spectrum of UHE protons 
can be approximated over a limited energy range by a power-law that fits 
the UHECR data below 60 EeV. This spectrum is then of the form A(z)E^ T 
where Ei is the energy of the proton. The UHECRs suffer energy losses from 
pair and pion production through interactions with the CBR and also cosmo- 
logical redshifting. The energy losses from pion production are determined 



according to equation (joj). The pair-production loss rate comes from [18 
In order to determine redshift loses, a flat ACDM universe with a Hubble 
constant of H = 70 km s _1 Mpc -1 is assumed, taking = 0.7 and Q m 
= 0.3. The source evolution is additionally assumed oc (1 + zp with Q = 
3.6, out to a maximum redshift of 2.5 which tracks the star formation rate. 
This value is a mean between the fast evolution and baseline models used in 



reference [19J. (See also references [20j, 12JJ.) The spectrum of UHECRs on 
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Figure 2: Comparison of the PAO data [9] with calculated spectra for various values of 
S^p. From top to bottom, the curves give the predicted spectra for S np = 1 x 10~ 22 ,6 x 
1CT 23 ,3 x l(T 23 , 1 x 1CT 23 ,0 (no Lorentz violation). 

Earth can then be determined from 

= 3cA(0) £ _ r /■*- + 

SvriJo Jo v/fi m (l + z) 3 + fi 

where A(0) is determined by fitting our final calculated spectrum to the 
observational UHECR data assuming T = 2.55, which is consistent with 
the spectrum derived by the Pierre Auger Observatory (PAO) collaboration 
below 60 EeVE The results of this calculation for various choices of the 
parameter 5^ p are shown in Figure [2] plotted along with the most recent 
results from PAO. 

It can be seen from Figure [2] that a small amount of LIV can still preserve 
the GZK suppression effect, but produces a "recovery" of the UHECR spec- 
trum at higher energies. Since this recovery is due to the virtual elimination 
of photomeson interactions at higher UHECR energies, it will also suppress 
the production of higher energy photomeson neutrinos. 



2 We have chosen a maximum UHECR energy of 3 x 10 21 eV. However, our results are 
insensitive to this value because they are determined by the LIV kinematics. 
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3. The Photomeson Neutrino Spectrum 



We now turn our attention to calculating the photomeson neutrino spec- 
trum that would result from the UHECR calculation detailed in the previous 
section. We use the data on the cross section for pion production compiled in 



reference [22[ and summarized in reference [23|. Near threshold the the total 
photomeson cross section is dominated by the emission of single pions. The 
most significant channel to consider involves the intermediate production of 



the A resonance 24 



p + 7 -»■ A -»■ iV + 7T (10) 

This channel strongly dominates the photomeson production process near 
threshold. Since the UHECR flux falls steeply with energy, it follows that 
the bulk of the pions leading to the production of neutrinos will be produced 
close to the threshold. 

For a proton interacting with the CBR, a pion and a nucleon are pro- 
duced. The outgoing nucleon has probability of 2/3 to be a proton and 1/3 
probability to be a neutron from isospin considerations. Should the resulting 
nucleon be a neutron, then the resulting pion is a tt + . Thus approximately 
twice the number of neutral pions are produced relative to charged pions from 
resonant pion production. However direct pion production, which accounts 
for about 20% of the total cross section, produces charged pions almost ex- 
clusively meaning that all told, approximately equal numbers of neutral and 
charged pions are produced around threshold. Neutral pions decay into pho- 
tons so we need only consider the charged pions for neutrino production. 
Three neutrinos of roughly equal energy result from the decay chain of the 

It is straightforward to determine the neutrinos produced and their ener- 
gies from the UHECRs. We follow closely the calculation of the neutrino flux 



as described in reference [12] and references therein. The key is to determine 
the amount of energy that is carried away by the pion. This follows directly 
from the inelasticity and the incident proton energy calculated using equa- 
tion OH]). We assume that all of the sources have the same primary injection 
spectrum and distribution as detailed in section [2j We calculate the total 
neutrino flux by integrating over proton energy, photon energy, and redshift, 
assuming the standard ACDM cosmology. We find that the shape of the neu- 
trino spectrum we obtain when assuming Lorentz invariance is very similar 



7 



-29 




Log E (eV) 



Figure 3: Neutrino fluxes (of all species) corresponding to the UHECR models considered 
in Figure[2] From left to right, the curves give the predicted fluxes for 5 np = 1 x 10~ 22 , 6 x 
1CT 23 ,3 x 1(T 23 ,1 x 1(T 23 ,0. 



to that obtained from the more detailed Monte Carlo calculations that in- 
clude all the relevant baryonic resonances and possible meson and multi-pion 
production channels (See, e.g., reference 13l|). 



The effect of LIV on the photomeson neutrino production is again man- 
ifested through the modification of the inelasticity of the interaction, since 
this determines the amount of energy that is carried away by the pion and 
therefore the resultant neutrino energy. The biggest impact of including LIV 
is to suppress the production of the higher energy photopions and therefore 
the resulting higher energy neutrinos. 

Figure [3] shows the corresponding total neutrino flux (all species) for the 
same choices of 5 np as the UHECR spectra presented in figure [21 As expected, 
increasing S np leads to a decreased flux of higher energy photomeson neutrinos 
as the interactions involving higher energy UHECRs are suppressed j3j . It is 
also evident that the peak energy of the neutrino energy flux spectrum (EFS), 
E<&(E), shifts to lower energies with increasing 5 np . We note that a similar 
effect can be produced on the photomeson neutrino spectrum from another 
possible mechanism involving LIV, namely "neutrino splitting" 25j]. With 



LIV, the decay of one neutrino into three neutrinos can be kinematically 
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allowed. This effect may also produce a decreased ffux at the high energy 
end of the neutrino spectrum. The feature that distinguishes between the 
two possible LIV effects is that neutrino splitting results in an increase in 
the flux of lower energy photomeson neutrinos. 



4. Considerations of UHECR Composition 



Throughout this paper, we have made the assumption that the highest 
energy cosmic rays, i.e., those above 100 EeV, are protons. The composi- 
tion of these primary particles is presently unknown. The highest energy 
events for which composition measurements have been attempted are in the 
range between 40 and 50 EeV, and the composition of these events is uncer- 



tain 



261, 27 , 28 



We note that in the case where the UHECRs with total energy above 
~100 EeV are not protons, both the photomeson threshold and the LIV 
effects are moved to higher energies because (i) the threshold is dependent 
on 7 oc E/A, where A is the atomic weight of the UHECR 29J, and (ii) 



it follows from equation (TTJ that the LIV effect depends on the individual 
nucleon momentum. We also note that the neutrino spectrum at the high 
energy end is the same for the mixed composition case as in the pure proton 



case 



31 31, 32 



5. Observational Prospects 

Several experiments currently place upper limits on the photomeson neu- 
trino flux. The ANITA long duration balloon experiment launched in De- 
cember of 2008 searched for electromagnetic cascades initiated by UHE neu- 
trinos within the Antarctic ice shelf via the Askaryan effect. Their analysis 
yielded a model-independent 90% CL limit on neutrino fluxes in the range 
of 10 18 - 10 23 eV with a sensitivity capable of excluding several optimistic 



cap 

photomeson neutrino flux models 15j. However, ANITA does not have 
sufficient sensitivity or energy range to distinguish a possible LIV effect on 
the neutrino spectrum. ANITA has an effective threshold energy ~ 10 18 
eV and its sensitivity about an order of magnitude too weak. IceCube is 
capable of detecting neutrinos < 10 17 eV, but its sensitivity is two orders 



of magnitude too weak to detect photomeson neutrinos [33]]. Proposed fu- 
ture space-borne missions such as the Extreme-Universe Space Observatory 
(EUSO) Q, super EUSO and Orbiting Wide-angle Light Detectors 
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(OWL) [36] would have much larger effective aperatures than presently avail- 
able detectors. They would be capable of making accurate determinations 
of the energy, arrival direction, and composition of the cosmic-rays and the 
associated photomeson neutrinos using a target volume far greater than is 
possible from ground-based experiments presently in operation. While such 
experiments could potentially provide the statistics necessary to observe the 
recovery of the UHECR spectrum at high energies, both EUSO and OWL 
as proposed would not achieve either the sensitivity or energies necessary to 
distinguish LIV effects in the neutrino spectrum. 

However, the proposed full Antarctic Ross Ice shelf ANtenna Neutrino Ar- 
ray (ARIANNA) would be capable of detecting photomeson neutrinos with 
a sensitivity an order of magnitude better than ANITA and other existing 
detectors. Like ANITA, ARIANNA exploits the Askaryan effect, i.e. the 
detection of coherent Cerenkov emission at radio wavelengths produced in 
the ice shelf by neutrino induced cascades. Because the power of coherent 
radio emission grows as the square of the shower energy and therefore neu- 
trino energy, ARIANNA is capable of detecting lower energy neutrinos than 
ANITA since the distances to balloon-borne detectors can be quite large 
while ARIANNA utilizes detectors situated directly on the ice shelf surface. 
ARIANNA is expected to observe ~40 events per 6 months in the en ergy 



range of photomeson neutrinos with energies in excess of ~ 10 17 eV 37 
This lower energy threshold is crucial for searching for LIV effects. As such, 
we shall restrict ourselves to the discussing the potential of ARIANNA for 
distinguishing the effect of LIV on the photomeson neutrino spectrum. We 
note that another proposed detector called IceRay would also make use of 
the Askaryan effect. IceRay would be placed at the location of IceCube and 



would also be capable of detecting photomeson neutrinos [38 



ARIANNA's high event rate, combined with its low energy threshold can 
distinguish LIV effects if 5 np < 3.0 x 10~ 23 with a 5 year exposure. We 
note that this limit is close to the upper limit indicated by the current PAO 
data Figure H] shows our calculated fluxes along with the sensitivity of 
ARIANNA for exposure times of 6 months 37J and 5 years. Here we have 
plotted E<&(E) for all neutrino flavors and v-v combinations compared with 
the proposed ARIANNA sensitivities to more clearly illustrate the threshold 
effect. We note that ARIANNA is sensitive to all neutrino flavors since it 
would primarily detect the hadronic shower and not the outgoing lepton. In 
Figure HI we present our neutrino fluxes and experimental sensitivities as the 
total of all flavors and v-v combinations. 
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Figure 4: All-flavor neutrino flux spectra (EFS) that correspond to the UHECR models 
considered in Figure [2] The six month sensitivity from the proposed AR1ANNA array 
is shown as a dashed curve. The dot-dashed curve shows the sensitivity scaled to a 
5 year exposure. From bottom to top, the solid curves give the predicted spectra for 
6n P = 1 x 1(T 22 , 6 x 1CT 23 , 3 x 1CT 23 , 1 x 10" 23 , 0. 
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It is clear from Figure H] that at the lower energy threshold of ~ 10 17 
eV for ARIANNA is very close to the peak energy in E<&(E) for 5 np = 3 x 
10 -23 and would have sufficient sensitivity at 10 17 eV to detect the expected 
neutrino flux provided it runs long enough to produce the desired sensitivity. 
ARIANNA is therefore very promising in terms of its ability to distinguish 
the effect of an amount of LIV that is compatible with the current PAO 
results. 



5.1. A Caveat 

We note that suggestions have been made to use the derivation of the ex- 
tragalactic 7-ray background flux up to 100 GeV using data from the Fermi 
Gamma-ray Space Telescope 



39 



to place limits on the absolute value of 
the photomeson neutrino flux, assuming that the 100 GeV 7- rays are pro- 
duced by an electromagnetic cascade off the CBR initiated by UHE pion 
decay 7-rays produced along with the neutrinos |40j, 41] . Such a constraint 



would lower the neutrino flux to almost an order of magnitude below the 
expected value and would make it much more difficult to test Lorentz in- 
variance with photomeson neutrinos. These results follow earlier neutrino 
constraints (42) obtained using various analyses of the EGRET-GRO results 
on the extragalactic 7-ray background 43], 44j. However, the argument here 
is contingent on the assumption that the extragalactic magnetic field is so 
small that a UHE electromagnetic pair-production-Compton cascade lead- 
ing to the production of 7-rays in the GeV energy range will not be cut off 
by synchrotron losses of the UHE electrons dominating over Compton losses 
At this point in time, the strength of the extragalactic B-field is 

9 G. 



45, 11 



46 . 47 



only constrained to be within the range ~3x 10 - ~ 3 x 10 

In addition, we note that the Fermi spectrum is not the result of direct 
observation, but of analysis. It critically involves the subtraction of both 
galactic foreground 7-rays and, in the case of Fermi instrumental calibration, 
by Monte Carlo modeling [48]. These are non- negligible uncertainties. As 
an example of the uncertainties involved, we note the significant differences 



431, 44] and the Fermi results [39(1 on 



between the EGRET-GRO results 
the extragalactic 7-ray background. 

We also note that since both the photomeson production cross section 
and the CBR photon spectrum are very well determined, and since the GZK 
cutoff effect is well documented [8j],[9j], a significant decrease in the predicted 
photomeson neutrino flux could either imply an unexpectedly small produc- 
tion of UHECRs or the existence of new physics. 
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6. Conclusion 

With future improved data from PAO, tighter constraints can be placed 
on the amount of LIV allowed by the UHECR spectrum. However even after 
a decade more of operation it seems unlikely that PAO would be able to deter- 
mine the UHECR spectrum with adequate statistics at energies greater that 
300 EeV, where the effect of LIV would be manifested. While space-borne 
missions such as JEM-EUSO and OWL could provide the statistics necessary 
to observe the effect, these missions are currently only in the planning stage 
and are many years from being realized. 

We have shown here that additional information on LIV or its constraints 
can be obtained from studying the spectrum of photomeson neutrinos. By 
calculating the flux of high energy neutrinos generated by the propagation of 
UHECR protons through the CBR in the presence of LIV, we find that LIV 
produces a reduction in the flux of the highest energy photomeson neutrinos 
and a reduction in the energy of the peak of the neutrino energy flux spectrum 
with both effects increasing with the strength of the LIV. Thus, observations 
of the UHE neutrino spectrum could provide a clear test for the existence 
and amount of LIV that would be exhibited in the highest energy cosmic-ray 
interactions. 

ARIANNA would have a sufficiently low threshold energy and the sen- 
sitivity necessary to determine the location of the energy peak in the pho- 
tomeson neutrino EFS to further test LIV in the range 8 np < 3 x 10~ 23 , 
consistent with the current limit indicated by the PAO data. The amount of 
LIV or its nondetection has important consequences for Planck scale physics 
and quantum gravity theories. 
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